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Abstract 

We have studied phase transition from hadron matter to quark 
matter in the presence of high magnetic fields incorporating the trapped 
electron neutrinos at finite temperatures. We have used the density 
dependent quark mass (DDQM) model for the quark phase while the 
hadron phase is treated in the frame-work of relativistic mean field the- 
ory. It is seen that the nuclear energy at phase transition decreases 
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with both magnetic field and temperature. A brief discussion of the 
effect of magnetic field in supernova explosions and proto-neutron star 
evolution is given. 
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1 Introduction 



The phase transition from hadron matter to quark matter (QM) 
is expected to occur in a variety of physically different situa- 
tions. In the astrophysical context this transition can take place 
in the interior of a cold, extremely dense neutron star. The 
other and the more likely situation in which such a transition 
can take place is the proto-neutron star (Prakash, Cooke, and 
Lattimer 1995; Prakash et al. 1997). These proto-neutron stars 
(PNS) have very high temperatures in their interiors (~50MeV) 
and more importantly, have a large amount of trapped neutri- 
nos (Burrows, and Lattimer 1986; Bethe 1990; Keil and Janka 
1995). It is well known that the structure of a PNS changes 
significantly in a few seconds after its formation as a result of 
the release of most of its energy and lepton excess by means of 
neutrino emission. After this PNS stage, the evolution of the 
neutron star (NS) may be described basically as a cooling pro- 
cess in which the thermal and structural properties are almost 
decoupled. It is worthwhile to look for signals of quark decon- 
finement during the short PNS evolution time scale, which is 
intimately related to the supernova event itself. Recently a lot 
of studies of the dynamical properties of core collapse on the ef- 
fect of the transition to QM have been undertaken (Mezzacappa 
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et al. 1998; Janka 1996). Such studies were motivated by the 
difficulties of the standard theoritical models in accounting for 
the observed explosions. In one of the recent studies Lugones 
and Benvenuto (1998) have considered the transition from nu- 
clear matter to QM under conditions prevailing in a PNS during 
the first tens of seconds since its birth. They considered a hy- 
peronic equation of state (EOS) for the nuclear phase and the 
'bag model' EOS for the QM phase. They studied the effect of 
the mass of the strange quark, m s , and the bag constant, B, on 
the phase transition. 

With the discovery of intense outbursts of low energy 7-rays 
from neutron stars there has been a lot of interest in the study 
of NS and their intense magnetic fields. Kuoveliotou et al (1998) 
have recently found a soft 7-ray repeater SGR 1806-20 with a 
period of 7.47s and a spin down rate of 2.6x10" syr^ 1 from 
which they estimate the pulsar's age to be about 1500 years and 
field strength of ~ 8xl0 14 G. Since most models of magnetic 
field generation give core fields to be 10 3 - 10 5 times higher than 
the surface fields, the magnetic field in the core of such neu- 
tron stars, called magnetars, could be as high as 10 18 - 10 19 G or 
even 10 20 G. Such intense magnetic fields may not be ruled out 
by the virial theorem if one takes significant general relativistic 
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corrections into account (Hong 1998). Further Kuoveliotou et al 
(1998) have estimated that upto 10% of all known NS could be 
these magnetars. Since it is well known that observationally it 
is very difficult to distinguish normal neutron stars from strange 
stars (SS) or hybrids (neutron stars with strange matter (SM) 
cores), some, if not all, of the observed intense magnetic field 
pulsars could be these exotic objects, formed by the deconfine- 
ment transition from nuclear matter to quark matter during the 
PNS stage itself under an intense magnetic field. 

Since such super strong magnetic fields affect the EOS (Anand 
and Singh 1999) and almost all known properties of NS or SS 
quite significantly, it will be interesting to study their effect on 
the deconfinement transition itself. Motivated by this, we in this 
paper study the viability of the transition from nuclear matter to 
quark matter in the conditions prevailing in such a magnetized 
PNS during the first few seconds of its birth. We follow closely 
the approach of Lugones and Benvenuto (1998, 1999). For the 
description of nuclear matter phase we have employed the well 
known mean field theory developed by Walecka, Serot, Glenden- 
ning and others (See Glendenning 1997 and references therein). 
We have included contribution of a, u, p interactions as well as 
scalar self interaction terms; with these one obtains the correct 
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experimental values of the various physical constants relevant 
in nuclear matter, viz, binding energy/nucleon, the saturation 
nuclear density, the asymmetry energy, the bulk modulus and 
the effective mass. Matter in this hadronic phase consists of 
neutrons, protons, electrons & electron neutrinos and their an- 
tiparticles; however we have not included the muon neutrino or 
the full baryon octet in this treatment. 

For the description of the QM phase, consisting of ud quarks, we 
have used the alternative description of QM in which the con- 
finement is treated by assuming a baryon density dependence of 
quark mass introduced by Fowler, Raha and Weiner(1981) and 
Plumber, Raha and Weiner(1984). Later on it was reformulated 
to show that the properties of SM in this density dependent 
quark mass model are quite similar to those predicted by the 
bag model and yet it is much easier to handle at finite tempera- 
tures. Recently we have used this DDQM model, in conjuction 
with strong magnetic fields, to study the radial pulsations of 
strange stars ( Anand et al. 2000). 

In section 2 a brief discussion of the EOS for the nuclear phase 
and the EOS for the QM at finite temperatures and in the pres- 
ence of a magnetic field is given. Also included in section 2 is 
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the criterion for finding phase transition from nuclear phase to 
QM phase described by Lugones and Benvenuto (1998). Section 
3 deals with results and discussion. 



2 The transition from nuclear matter to quark 
matter 

2.1 The equation of state for nuclear matter 

Following Lugonese and Benvenuto (1998), we apply the nuclear 
mean field model for nuclear matter EOS. We incorporate n, p, 
e and v e particles and their antiparticles (since the antiparticles 
are also present at a finite temperature). In the absense of mag- 
netic field, expressions for the number, pressure and energies 
densities are 

n ' = ~(2^f d3p[MT) ~ MT)] (1) 

1 g, , p 2 d 3 p[fj(T) + j,(T)\ m 
' 3(271-)^ (pt + mlyl 2 y > 

pi = WfI £p{p2 + m * 2 ) 1/2 [^ T ) + ^ ( 3 > 

where i=n, p, e, v e . The spin degenracy factor, gi, is unity 
for v e and two for all others. Here fi(T) and fi(T) are the 
fermi-dirac distribution functions for particles and antiparticles 
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respectively: 

fi(T) = (exp([^ 2 + m2 - ^]/T) + l)" 1 (4) 



/i(T) = (exp([vWm? + ^]/T) + l)" 1 (5) 

These expressions are valid for all the particles; however for 
neutron and proton, as a result of the strong interactions as de- 
scribed by the mean field theory, the mass m 8 has to be replaced 
by the effective mass m* and the chemical potential pi by the 
effective chemical potential p*: 

m* = m l - g a a (6) 



2 ^ 

Ph =IH- (— J - — P/ 3 (7) 
Vrrw \m p J 3 

The iso-spin density p lz and the baryon densities p B are given 

by 

Ph = l^pp-Pn) (8) 
Pb = Pp + Pn (9) 

The total pressure P and the mass energy density p for the 
system are given by 

- v<** + k(£f * (10) 



+ \c(g a aY + U^\ pi (11) 



I w ' 2\m Pl 



where (Glendenning, 1997) 

(^X = 11.79 fm~ 2 , (^X = 7.U9fm- 2 
\m a J Km^J 

\ 2 

^ =4.411/m" 2 , 6 = 0.002947, c= -0.001070(12) 

The mean field g a a satisfies 

(—) {9* <j) + bm N (g a a) 2 + c {g a af = En? (13) 
\m a J 

Here g a , g^, g p , m a , m w , m p are the sigma, omega, rho couplings 

and masses respectively, b and c are the strengths of the two 

scalar self interaction terms in the Lagrangian: 

— % 
L = E ipi[i>7^ + + -g p r .p - {m-g a a)]ipi 

i=n,p ^ 

+\(d,ad»a -ml a 2 )- + \mlu^ 

In addition to the number density for neutron and proton we 
also require the scalar density n\: 



9i f rn*d*plfi(T) + /;(T)] 



ttW (v 2 + mfY/ 2 [ } 



(2nfJ (p 2 + mf)V2 
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We now introduce a constant and uniform magnetic field of mag- 
nitude B in the positive Z direction. The motion of a charged 
particle of charge e and mass m in such a magnetic field gets 
quantized and the energy eigenvalues are given by 
E 2 = p\ + m 2 + 2(n + 1/2 + s)eB 

where n=0, 1, 2, 3, defines the various Landau levels and 

s=±l/2 refers to the spin up (down) states of the particle. The 
number of states per unit volume in an interval of longitudinal 
momentum Ap z for a given level n is (2 — 6 no )-^0^Ap z . [we are 
using h = c = 1] . Hence the usual sum over particle states (per 
unit volume) in zero field : 

(2 being the spin degenracy factor) must be replaced by 

^E n (2-Sno)fd Pz 

With this ansatz it is straightforward to obtain expressions for 
the energy density p, pressure density P, baryon density n and 
the scalar density n s of a charged particle in magnetic field. 
These expressions are: 

ni = ^ E(2 " « / dp z [fi(T) - UT)] (16) 

a 9ieBm*i r dp z [fj(T) + fj(T)] 
n, = E(2 " Sno) j ^ (17) 



Pl = ^V 2 ~ 6no) J Ej (18) 
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Pi = 



2tt 2 



E( 



2 - / Eidp z (fi(T) + /i(T)) (19) 



The hadron phase is assumed to be charge neutral and in j3- 
equilibrium. Electric charge neutrality requires 



Chemical /^-equilibrium in the presence of trapped electron neu- 
trinos implies the following relation between the chemical po- 
tentials of various particles involved in the hadron phase: 



All the above equations can be solved numerically by giving four 
quantities, the temperature T, the electron chemical potential 
/i e , the chemical potentials of the electron neutrinos ji Ve and the 
magnetic field B. 

2.2 The equation of state for quark matter 

The quark phase is composed of u, d quarks, gluons, electrons 
and electron neutrinos. We shall treat QM as a free fermi gas 
of u, d quarks where the mass of each quark is parametrized as 



Tip — TL e 



(20) 



(21) 



C 



(22) 



m u = m d = 



3n B 
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C being a constant. The range of C is to be constrained by 
stability arguments. The pressure in the quark phase is given 

by 

P = EP 4 (23) 
where the sum goes over i = u, d, e, v e and their antiparticles. 
For the electron and the electron-nentrino, pressure is given by 
eq(2) in the absence of magnetic field and by eq(18) in the pres- 
ence of magnetic field. To obtain expressions for P, p and n for 
u and d quarks, we use the thermodynamic potential which in 

the presence of magnetic field, B, is given by 
ap-BT r 

^ = -%r- ?( 2 " <W / <fr*[iog(i + eM-m - *])) 

+ log(l + exp(-/?[£? + ^]))] (24) 
where 



E* = yffi + mf + 2nq % B (25) 

gi= 3 (colour) and (2/3,-1/3) for u and d quark respectively. 

The pressure Pi is given by (Peng et al. 2000) 

dQ 

Pi = n B —^ - Qi (26) 
on B 



energy density pi by 



dQi 

Pi = Qi + piUi - T— (27) 



and number density ni by 



dQi , x 

n % = 28 

dm 
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Using eq(24) for the thermodynamic potential f^, yields: 

SieiB„,„ r , r pldp 2 [fi(T) + fi(T)} 



2tt 2 V v '""•/ 

C 9i eiBmi^ fn . f dp z [fi(T) + /,(T)] 



E(2-U/ ™^^ (29) 



3n B 2tt 2 V v ^ £?? 

The expressions for number density and energy density remain 
the same as in hadron phase. 



2.3 Transition from PNS to quark matter 



Following Benvenuto et al it is assumed that the transition from 
nuclear matter to quark matter in chemical equilibrium occurs 
necessarily through an intermediate step of quark deconfinement 
driven by strong interactions. In view of this, in a deconfinement 
transition the abundance per baryon of each quark and lepton 
flavour is the same in both the hadron and quark phases. This 
is in contrast to the model first given by Glendenning (1992) 
who proposed the existence of a mixed phase. 



In order to compute the conditions for phase transition, we 
apply the standard Gibb's criteria i.e, equality of pressure P, 
temperature T and Gibb's energy per baryon g in both phases: 

P q = Ph, T q = T h , g q = g h (30) 
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together with flavour per baryon coservation equations 

= y? (3i) 

for i=u, d, e, v e . Notice that this condition automatically makes 
the quark phase to be charge neutral. Here 

yx _ J}i_ 

n B 

for x = h or q. Also 

Yu = ^Y p + Y r , 
Y d h = Y p + 2Y r , 

Starting from Y p and Y n , obtained from the hadronic phase, one 
then obtains various number densities and thereby the chemi- 
cal potentials in the quark phase; energy and pressure are then 
obtained from Eqs(19) and (29) respectively. 



(32) 



(33) 



3 Results and discussion 

In this work we are interested in the deconfinement transition 
from hadron to quark matter in the presence of trapped neu- 
trinos and magnetic field. As discussed in detail by Benvenuto 
and Lugonese, we see that for a fixed set of values of n„ e , T 
and eB the nuclear phase depends only on one of the chemi- 
cal potentials (say /i p ). With /j, p , /i e , T and eB fixed, we can 
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calculate n p , n n , n e , n Ve and so the abundances Y p , Y n , Y e , Y Ve , 
Y*, Y£. Using conditions (30), (31) we obtain the mass energy 
density U of hadron matter at which it deconfines. At finite tem- 
perature, we have to evaluate numerically Fermi-Dirac integrals 
appearing in nf, Pi and E{. We were able to evaluate these 
integrals upto an accuracy of about 0.1% by using a technique 
due to Cloutman (1989). At this stage it is necessary to explain 
this method by means of a test integral. Consider the following 
integral occuring in nf. 

J—oo 

= 2 r % ., (34) 

Here 

^* = (rf + < 2 ) 1/2 (35) 

Putting j3p z =x, we have 

roo dx 
Ii = 2T - , 2 , 2u/2 — T (36) 

J0 14- e (x 2 +m i 2 ) 1 / 2 -iJ, i v J 

where 

m'i = m*(3 

Hi = tip (37) 

We can write another form for Ii by substituting 

t= (x 2 + m' 2 ) 1/2 - m'i (38) 
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VIZ., 



L = 2T 



where 



oo (t + m\)dt 



yt(t + 2m' i )][l + et-v 



= 2 T r-^)dt (39) 



4>if) = , t + mi , (40) 

't(t + 2m-) V J 



Depending upon the values of r] we can have various forms of I{ 



Case 1: rj > 30 

J< = m + 2 £ C 2j <f> 2j (41) 

where 

^ - (42) 

Case 2: -5 < r/ < 30 

The original form (eq 34) is good for evaluating Ii numerically 
by Simpson's rule. The upper limit and the number of points 
were adjusted in such a way that the evaluation of integral was 
accurate to about 0.1%. 
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Case 3: r\ < — 5 

We expand U in terms of modified bessel function K\. 

li-H (-l)Xe M>+1) ^iK(^ + 1)) (43) 

Using the above method of evaluating integrals at finite tem- 
peratures, the transition from PNS to QM can be determined 
at a given temperature for various values of fi^ provided the 
parameter C is known. We choose C in such a way that the ud 
system is unstable. We have chosen three sets of values of C 
namely 90, 110, 130 MeV/m" 3 . 

In fig 1, we have plotted the mass energy density, U, of hadron 
matter (in terms of the nuclear saturation density, p = 2.7 
x 10 14 gmcra~ 3 ) at which phase transition occurs verses the chem- 
ical potential of neutrino present in hadron matter for three val- 
ues of C, viz, C=90, 110, 130 MeV/m -3 and at temperatures 
T=5 and 60 MeV. We find that the phase transition density is 
an increasing function of the neutrino chemical potential. It is 
interesting to note that the phase transition density increases 
with increase in the value of C. On the other hand, the energy 
density is a decreasing function of temperature T. Our results in 
this regard are similar to those of Lugones and Benvenuto (1998) 
who have used the bag model EOS for the quark system. In both 
cases the energy density required for the transition to take place 
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increases with increasing p, v ^ but for given p, Ve decreases with 
increasing temperature. In quantitative terms the energy den- 
sity at transition in our model is less by about a factor of two 
as compared to the Lugones-Benvenuto results. In fig2 we have 
studied the phase transition of hadron matter to quark matter 
in the presence of a strong magnetic field. We have taken mag- 
netic fields of 0, 5xl0 4 and 10 5 MeV 2 (1 MeV 2 ~ 1.6xl0 14 G) at 
temperatures 5 and 60 MeV and C=90 MeV/m~ 3 . The results 
are shown as a graph of U vs. p, Ve . It is interesting to note that 
the phase transition density decreases with increase of magnetic 
field at a given temperature. Similar results are obtained for 
higher C as well at a given temperature. Figure 3 deals with 
the behaviour of phase transition density with temperature for 
different values of p Ve . The phase transition density decreases 
from 2.2p to about l.Sp Q for eB=0 when \i Ve varies from 300 
MeV to MeV, whereas for eB=10 5 MeV 2 U falls from about 
2p to 0.8/v Figure 4 shows the same kind of results for C=110 
MeV/m- 3 . 

To conclude, we find the phase transition density is lowered in 
the presence of magnetic field. In this regard the effect of mag- 
netic field is similar to that of temperature. The presence of a 
strong magnetic field is likely to cause a transition to free quark 
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matter in the proto-neutron stars in early stages of its evolution. 
The question can however be settled by studying the EOS and 
M-R relation for such a magnetized PNS to see whether it meets 
the conditions of the energy density required for the transition 
to take place. This question is under study and will be reported 
elsewhere. 
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Figure captions 

Figure 1. The mass energy density, U, of hadron matter at 
which phase transition occurs versus the chemical potential of 
the electron neutrinos present at eB=0. Labels al, a2 and a3 
correspond to T=5 MeV and C=90, 110 and 130 MeV/m" 3 
respectively, whereas bl, b2 and b3 correspond to T=60 MeV 
and C=90, 110 and 130 MeV/m" 3 respectively 

Figure 2. The mass energy density, U, of hadron matter at 
which phase transition occurs versus the chemical potential of 
the electron neutrinos present at C=90 MeV/m -3 . Labels al, 
a2 correspond to eB =0 and T=5 and 60 MeV respectively; 
bl, b2 correspond to eB=5xl0 4 MeF 2 (1 MeV 2 ~ 1.6xl0 14 G), 
T=5 and 60 MeV respectively and cl, c2 correspond to eB= 
10 5 MeV 2 T=5 and 60 MeV respectively 

Figure 3. The hadron matter mass energy density of phase 
transition, U, versus the temperature T at C=90 MeV/m~ 3 . La- 
bels aO, al, a2, a3 correspond to eB=0 and ^^ e =0, 100, 200, 300 
MeV respectively; bO, bl, b2, b3 correspond to eB=10 5 MeT /2 
and fj, v =0, 100, 200, 300 MeV respectively 

Figure 4. The same as fig. 3 but for C=110 MeV/m~ 3 . 
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Figure 1: The mass energy density, U, of hadron matter at which phase tran- 
sition occurs versus the chemical potential of the electron neutrinos present 
at eB=0. Labels al, a2 and a3 correspond to T=5 MeV and C=90, 110 and 
130 MeV/m~ 3 respectively, whereas bl, b2 and b3 correspond to T=60 MeV 
and C=90, 110 and 130 MeV/m^ 3 respectively. 
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Figure 2: The mass energy density, U, of hadron matter at which phase tran- 
sition occurs versus the chemical potential of the electron neutrinos present at 
C=90 MeV/m -3 . Labels al, a2 correspond to eB =0 and T=5 and 60 MeV 
respectively; bl, b2 correspond to eB=5x 10 4 MeV 2 (1 MeV 2 ~ 1.6xl0 14 G), 
T=5 and 60 MeV respectively and cl, c2 correspond to eB= 10 5 MeV 2 T=5 
and 60 MeV respectively. 
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Figure 3: The hadron matter mass energy density of phase transition, U, 
versus the temperature T at C=90 MeV/m~ 3 . Labels aO, al, a2, a3 corre- 
spond to eB=0 and /v=0, 100, 200, 300 MeV respectively; bO, bl, b2, b3 
correspond to eB=10 5 MeV 2 and fi Ve =0, 100, 200, 300 MeV respectively 
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